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ABSTRACT 

The maximum efficiency of silicon solar 
cells is calcifLated and yields a value near 18 $. 
Additionally, the performance of these high effi- 
ciency cells in a synchronous orbit radiation 
field is calculated and indicates that these 
cells would be superior to present silicon cells 
at all times. The performance of conventional 
cells is analyzed and several areas in which 
performance gains may be achieved are discussed. 

These, areas .include improvements in diffused 

region profile, in reduction of excess "forward 
currents in cells made from low resistivity ma- 
terial and in the theory for describing complex 
solar cell structures. 

INTRODUCTION 

Over the last several years, interest in the 
maximum efficiency of silicon solar cells has 
been reawakened. Previous studies (1, 2) con- 
cluded a maximum efficiency of 22% and indicated 
several areas in need of study before this effi- 
ciency could be achieved. It is the purpose of 
this paper to repeat these efficiency calcula- 
tions with more realistic loss terms due to re- 
flection and grid coverage and a more general 
expression for the reverse saturation current. 
Next, because cells are widely used in a space 
environment, the performance of these improved 
cells in a synchronous orbit radiation field was 
calculated. Finally, the performance of present 
cells was analyzed in terms of present theories 
to determine areas which hold the most promise 
of yielding significant gains in efficiency. 


decreasing recombination at the surface of the 
cell, and by decreasing both surface reflection 
and front contact shadowing. The fill factor 
can be raised by reducing the series resistance 
and improving the junction characteristics. In 
the ensuing pages attention will be directed to- 
ward a*l 1 these avenues to raising efficiency. 

Efficiency Calculations 

Diode Equation: The starting point for the 

ultimate efficiency calculations was the simple 
— diode equation: •_ 

I = I 0 (e kT -l)-I L (1) 


where I is the current, I Q is the reverse satura- 
tion current calculated from diffusion theory, 

V is voltage, q is electronic charge, k is 
Boltzmann's constant, T is the absolute tempera! 
ture (°jc) and 1^ is the light generated current. 
Because diffusion theory is being used, the A 
value generally shown in equation (l) equals 
unity and therefore can be eliminated. 


Open-Circuit Voltage: I Q is the parameter 

in equation (l) which controls the open- circuit 
voltage. Therefore, its magnitude and variation 
with material parameters is important „ The gen- 
eral form of for any base width is as follows: 
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THEORY AND CALCULATIONS 

The purpose of this section is to determine 
the maximum practical efficiency of a silicon 
cell. The power output, F max of the cell (and 
hence its efficiency) can be written in the fol- 
lowing way: 


where F.F. is fill factor, I sc is the short-cir- 
cuit current and V QC is the open-circuit voltage. 
Raising any of the three quantities increases the 
power output. In general, the open- circuit volt- 
age can be increased by reducing the resistivity. 
The short-circuit current can be increased by in- 
creasing the base minority carrier lifetime, by 


where n. is the intrinsic carrier concentration 
in silicon, N^Np) is the acceptor (donor) con- 
centration, D n (D_) and L^Lp) are the diffusion 
constant and the^diffusion length of electrons 
(holes) in the p-(n-) region, W (W n ) the width of 
the base region (diffused region) of the device. 
The assumptions inherent in equation (2) include 
an infinite surface recombination velocity at 
each face and uniform doping density in each re- 
gion. These latter two assumptions are known to 
be incorrect for conventional solar cells, and at 
the present time, it is not possible to estimate 
the magnitude of error that they introduce. 

Materials constants were obtained from the 
literature (1-7). The cell thickness was fixed 
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at 300 pM and the -width of the diffused region at 
0.25 uM. Base minority carrier lifetimes were 
obtained from lies (2) for resistivities greater 
than 0.01 fi-cm, Wolf's data were used for 0.01 Or- 
em. The lifetime in the diffused region was 10”® 
sec. Table I summarizes the data used and the 
results of the calculation of Iq. The parameters 
chosen for the diffused region (i.e. Kp = 5 x lO^, 
tL = 20 cm 2 /V-sec, t = 10"° sec, W = .25 uM) yield 
a p fixed value for the I Q from the diffused region 
of 1.2 x 10" " A/cm 2 . This leads to a value of 
I Q in 0.01 fi-cm material higher than calculated 
in references 1 and 2. 

Light Generated Current: Wolf (8) has de- 

rived a general equation for minority carrier 
current collected at a n-p junction under illumi- 
nation. • This equation has been programmed for 
use on a digital computer. A copy of the program 
was obtained from Wolf and modified slightly to 
allow additional flexibility. These modifica- 
tions ijiclude calculation of depletion region 
width using an 'abrupt junction approximation -and - - 
the inclusion of spectral reflectance of the sur- 
face so that spectral response curves can be cal- 
culated as part of the output data. 

Values of the light generated current were 
calculated from this program as a function of 
both base resistivity and base diffusion length 
for a standard n + -p configuration. The surface re- . 
combination velocity was chosen to be 10 2 cm/sec 
at the front of the cell and 10° cm/sec at the 
ohmic rear contact. A complementary error func- 
tion profile was assumed in the diffused region. 

In the program, this impurity profile was approx- 
imated by two exponential fields. Base doping 
was uniform and cell thickness was 300 pM. In- 
clusion of a p + drift field at the back contact 
had only a small effect (~ 2%) on I for 10 and 
1 fi-cm resistivities. However, because the equa- 
tion for 3^ does not include the effects of drift 
fields, they are not presented in this paper. 

All calculations were performed assuming a total- 
ly absorbing silicon surface and the resulting 
currents were reduced by 5 % for a hypothetical 
uniform antireflecting surface. 

Figure 1 shows the variation of short-cir- 
cuit current with diffusion length calculated 
with the parameters listed above. There is lit- 
tle increase in current above diffusion lengths 
of 200 uM as has been noted before in silicon 
cells. For a constant diffusion length, the cur- 
rent decreases as resistivity decreases because 
the depletion region gets smaller. 

Series Resistance: Series resistance is un- 

avoidable in the conventional solar cell, largely 
due to the sheet resistance of the diffused layer 
and the grid finger spacing ( 9 ). Measurements 
of series resistance on cells using standard 
techniques (10, 11) yielded series resistances 


about 0.3 0. This resistance will lower the 
power by 3%. 

Grid coverage reduces the active area of a 
solar cell and hence its power output. Present 
cells have grid coverages between 5 and 10%. The 
lowest area coverage would came from a wrap- 
around geometry. Using standard 6-line coverage 
on a 2 x 2 cm cell the wrap-around geometry would 
lead to a covered area of 5%. For these calcu- 
lations then, the cell current was reduced by 5% 
to account for contact coverage. 

Efficiency Calculations: Equation (1) can 

be solved for short-circuit current, I EC , open- 
circuit voltage, Voc > and the maximum power 
point, p . The solutions for I sc and V oc are 
straight?orward. The solution for maximum power 
results in a transcendental equation that can be 
easily solved by iteration using the Newton- 
Raphson technique. This solution is especially 
suited for modern desk-top calculators with mem- 
ory. - -Value S- of -I Q calculated from equation (2) 
and values of 1^ obtained from figure 1 were' 
used. The values of I L on figure 1 were re- 
duced by 5/ to account for the coverage of the 
grid as noted above. The re suiting values of 
and efficiency were then reduced by 3 % to 
account for a 0.3 £1 series resistance. The val- 
ues of V Q _, I EC and efficiency so calculated are 
shown in figure 2. The peak efficiency is about 
18$ occurring for a base resistivity between 0.1 
and 0.01 fi-cm. Table II lists the performance 
obtained for the various resistivities. The 
maximum efficiency is slightly lower than that 
reported by Wolf (1) because of higher loss val- 
ues chosen for reflectivity ( 5 % to his 3%) and 
grid area coverage (5$ to 0%) and the higher I 
value calculated for 0.1 fi-cm material. 

In summary, for a simple planar p-n junction 
model, a practical upper limit of silicon solar 
cell efficiency of 18% has been obtained for 
0.1 fi-cm material. 

Radiation Damage 

Diffusion Length Degradation: Estimates of 

the perfomance of these high efficiency cells in 
a radiation field were made using the standard 
model for diffusion length degradation: 



where L is the diffusion length after same flu- 
ence $, L is the initial diffusion length and 
K is a damage coefficient for the bombarding 
particle and the semiconductor material. For 
this calculation synchronous orbit was assumed 
as was the damage integral (K $) as calculated 
by Cooley and Barrett (12) for 10 fi-cm material _ 
covered with .015 cm fused silica (K$ = 4 . 12 x 107 
cm 2 - year). It was assumed that the cover glass 
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used, did not lower the initial short-circuit cur- 
rent. The damage integral so calculated assumes 
the damage to be caused by solar flare protons 
and trapped electrons. 


The problem of determining damage integrals 
at other resistivity values is a bit more diffi- 
cult. From data of Denney and Downing (12, 13, 
14), it was estimated that the damage integral 
for 1 fl-cm silicon material was about 3-5 times 
that of 10 fl-cm material. For this calculation^ 
the worst case was chosen (i.e., KJ> = 2.06 x 1 Cr/ 
cm 2 -yr). Below 1 fl-cm material, however, almost 
no data are available. Curtis and Srour (15) pre- 
sent data on neutron-irradiated material which 
indicate that the damage coefficient becomes con- 
stant from 1 fl-cm to lower resistivities. Be- 
cause there is a general correspondence between 
neutron and proton or electron damage in silicon, 
the damage integral was assumed to remain constant 
at 2.06 x lO^/cm -yr for resistivities 1 fl-cm and 
below. 


Performance Calculation: The change in dif- 

fusion length determined from equation (3) was 
used to calculate both an I from equation (2) 
and a light generated current from figure 1. These 
data were then used to determine cell parameters 
using equation (1) and the corrections previously 
noted. Figure 3 shows the results of these calcu- 
lations. Also included are the calculated values 
for a present-day 10 fl-cm cell and same data ob- 
tained after 1.1 and 3.3 years ( 1 6 , 17) in syn- 
chronous orbit for conventional cells. These lat- 
ter data demonstrate the general validity of the 
calculations. It can be seen that a.n of the 
"advanced" cells are significantly higher in effi- 
ciency than the present day cell at all times. 

The advantage comes largely from better collection 
efficiency, optical coatings and fill factor. Due 
to its higher damage coefficient, the efficiency 
of the 1 fl-cm cell falls below that of the 10 fl-cm 
cell after about 5 years in orbit. 

PATHS TO IMPROVE EFFICIENCY 

Characteristics of present solar cells were 
examined in the light of the same theory and as- 
sumptions used to calculate the ultimate effi- 
ciency. The purpose was to assess what areas of 
improvement offer the most gain in efficiency and 
what knowledge must be acquired to realize improve- 
ments in these areas. 

Improvements in Short-Circuit Current 

As can be seen from figure 1, increasing the 
diffusion length from 150 to 600 pM in the silicon 
solar cell produces only about a 1 mA/cm? increase 
in cell output. Therefore, greatly increased dif- 
fusion lengths will not greatly increase cell out- 
put. Furthermore, the diffusion lengths in 0.1 fl- 
cm material are on the order of 90 pM (t = 4 psec) 
(18), leading to currents within 10$ of the maxi- 
mum achievable. Hence, major efforts to increase 


the diffusion length in this material would not 
produce large increases in output. 

The measured spectral response of a 10 fl-cm, 
uncoated, 0.25 pM junction silicon solar cell is 
shown in figure 4. An insight into same of the 
materials parameters can be obtained by fitting 
this response with the theory described previous- 
ly. The response of the cell in the infrared 
spectral region is strongly dependent on diffu- 
sion length. A value of 170 pM was obtained on 
the data shown. Direct measurement of the diffu- 
sion length by the X-ray technique (19) yielded 
a value within 5$ of that obtained by fitting. 

As can be seen, however, it was not possible to 
fit the blue response of the cell by assuming a 
complementary error function for the impurity 
profile in the diffused region. Only by using 
the anomalous profile noted by lies and Liebenhaut 
(20) could satisfactory agreement be obtained. A 
front surface recombination velocity of lo5 cm/sec 
was used. Subsequent measurements at NASA-Lewis 
_of the impurity profiles in similarly processed 
cells have indicated the presence - of an anomalous - 
profile . 

The response of a cell to blue wavelengths of 
light is significantly reduced when both an anom- 
alous profile and a high surface recombination 
velocity are present. An example of this effect 
is shown in figure 5, where the calculated re- 
sponse at 0.4 pM of a SiO-coated silicon cell is 
shown for various values of S and diffused im- 
purity profile. It appears that a significant 
increase in output of a solar cell could he a- 
chieved by obtaining a complementary error func- 
tion profile even if S were to remain high. Cal- 
culations show that this increase would amount to 
nearly 3 mA/cm 2 . Conceivably, if an error func- 
tion profile were achieved, S might also he re- 
duced because of the resulting improvement in 
perfection of the surface; however, subsequent 
reduction of S to 10 2 cm/sec would add less than 
0.5 mA/cm 2 to the cell output. 

Improvements in Open-Circuit Voltage 

The key to obtaining significant increases 
in solar cell efficiency lies in obtaining large 
increases in the open-circuit voltage by reducing 
the resistivity. The open-circuit voltage is 
strongly related to I Q (for a given I Tj ) . Simple 
diffusion theory was used to predict the behavior 
of cell open- circuit voltage for a conventional 
10 fl-cm cell with temperature. The results of 
this calculation are shown in figure 6. Perform- 
ance data were taken from previously published 
work (21). Since excellent agreement is obtained, 
it appears that the open-circuit voltage of 10 fl- 
cm cells can he adequately described by simple 
diffusion theory with unity A values. 

It has, however, been noted by lies (2), 
that the open-circuit voltage of cells does not 
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continue to increase as the resistivity decreases 
as expected frcan diffusion theory. Instead, the 
open-circuit voltage saturates at about 0.62 V 
for resistivities lower than 1 fl-cm. This volt- 
age difference cannot be attributed to poor bulk 
minority carrier lifetimes. lies (18) has found 
that 0.1 ft- cm cells with diffusion lengths of 
90 mM can show open-circuit voltages in the 0.55- 

0.6 V range. To better understand the reason for 
this loss in voltage, the forward diode character- 
istics were obtained on both 10 and 0.1 fl-cm 
cells. lies also etched small mesa regions (10 
to 20 mesas each with a diameter of 0.145 cm per 
wafer) into several of these 0.1 and 10 ft-cm 
cells to sample the distribution of characteris- 
tics over the surface. The results of these meas- 
urements are shown in figure 7 where envelopes 
encompassing the typical range in performance for 
both resistivity materials are shown. A line at 
a current density of 35 mA/cm 2 (equivalent to 
sunlight illumination) is provided for reference. 
The large excess .forward currents noted at low 
forward voltages are high enough to significantly 
reduce the open-circuit voltage in 0.1 fl-cm ma- 
terial below that expected. The variation from 
spot to spot is also quite large. In 10 ft- cm ma- 
terial both the magnitude and variability of the 
excess forward current is less than for the 0.1 Q- 
cm case. Accordingly, there is no loss in open- 
circuit voltage over that predicted. (There is a 
small loss in maximum power, however. ) 

The mechanism giving rise to this excess for- 
ward current is not clear at the present time. 
Possible explanations for this effect include cur- 
rent loss through defect regions in the diffused 
layer and junction and surface or edge leakage 
currents around the junction. The increase in 
efficiency predicted earlier depends strongly on 
increased open- circuit voltage in low resistivity 
material. Therefore, it is imperative that the 
cause of this excess forward current be under- 
stood and eliminated so these gains can be achiev- 
ed. 

Improvements in Theory 

Increased open-circuit voltages (.58-, 585 V) 
have been reported for n + -p-p + solar cell struc- 
tures made from 10 ft-cm material (2, 22, 23). 

These voltages are significantly in excess of the 
0.55 V maximum open-circuit voltages seen in con- 
ventional cells. The use of simple diffusion 
theory to obtain I Q as outlined previously does 
not explain these high voltages, even when long 
minority carrier lifetimes are used. Equation (2) 
was derived on the basis of uniform doping in 
both regions. Therefore it is not entirely sur- 
prising that it does not predict the voltage in 
the n + -p-p + case, where significant doping varia- 
tions are present. Eefore these higher open-cir- 
cuit voltages can be explained, the diffusion ex- 
pression for I must be rederived for the more 
general case with arbitrary doping profiles. The 
further insight gained from this work may lead to 
further increases in cell output. 


COMMENTS AND CONCLUSIONS 

Calculations indicate that a potential in- 
crease in efficiency of silicon solar cells to 
near 18% may be achievable through the use of 
low resistivity material possessing good diffu- 
sion lengths, diode characteristics and improved 
surfaces and diffusion profiles. Further, the 
calculations indicate that these cells should be 
superior to present day silicon cells in a syn- 
chronous radiation environment. 

Analysis of conventional cells, however, dis- 
closes a number of problem areas that must be 
studied further in order to achieve higher out- 
puts. The increases in short-circuit current 
that may be achievable in the cell are rather 
modest. Because the diffusion lengths in present 
cells (both 10 and 0.1 ft-cm) are long, only minor 
increases in current could be expected from large 
increases in diffusion length. Increases in the 
blue response of cells can be achieved by elimina>-. 
tion of an anomalous impurity diffusion profile 
in the diffused region. If a complementary er- 
ror function profile were obtained, calculations 
show that an increase of nearly 3 mA/cm 2 could 
be achieved even if the surface recombination 
velocity were to remain high. 

The measured open- circuit voltages of low 
resistivity material are much lower than pre- 
dicted by diffusion theory. This loss is not 
due to poor diffusion lengths. Rather it is 
due to excess forward diode currents. In order 
to achieve the large anticipated increase in 
open- circuit voltage for low resistivity mater- 
ial, the origin of this excess current must be 
determined and eliminated. 

Although simple diffusion theory yields open- 
circuit voltage in agreement with practice for 
conventional 10 ft-cm n+-p cells it does not yield 
satisfactory results for more complex structures 
such as the n + -p-p + cell. Increased theoretical 
understanding of these cells may lead to further 
insights that may yield additional improvements 
in the output of solar cells. 

Thus a variety of problem areas that limit 
increased output from silicon solar cells has 
been outlined. An agressive program aimed at 
these problems can be expected to either yield 
increased output or at least an understanding of 
the ultimate limits of efficiency in the silicon 
solar cell. 
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CALCULATED SHORT CIRCUIT CURRENT DENSITY 
AS FUNCTION OF RESISTIVITY 
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Fig. 1. 


CALCULATED PERFORMANCE OF SILICON 
SOLAR CELLS AS FUNCTION OF RESISTIVITY 



RESISTIVITY, O-CM cs-am 

Fig. 2. 
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CALCULATED DAMAGE TO ADVANCED SOLAR CELLS 
IN SYNCHRONOUS ORBIT 


VARIATION OF OPEN-CIRCUIT VOLTAGE OF 
lOfl-CM N/R SILICON SOLAR CELL 
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SPECTRAL RESPONSE, OF lOfl-CM 
SILICON SOLAR CELL 


FORWARD DIODE CHARACTERISTICS OF MESAS 
MADE ON 10 AND O.lfl-CM 
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CALCULATED SPECTRAL RESPONSE AT 0.4/jM OF 
SiO COATED SILICON SOLAR CELL 



Fig. 5. 
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TABLE I 


CALCULATION OF SATURATION CURRENT DENSITY FOR n + -p SILICON SOLAR CELLS 


RESISTIVITY 

MAJORITY CARRIER 
CONCENTRATION 
N. cm -3 

A, 

MINORITY CARRIER 
MOBILITY 
p, cnr/V-sec 

MINORITY CARRIER 
LIFETIME (BASE) 
T 5 sec 

SATURATION 
CURRENT 
DENSITY 
V A /cm 2 

10 

1.35 x 10 15 

1300 

-4 

1 x 10 

2.64 x lo' 11 

1 

1.6 x io l6 

1100 

8 x lo" 5 

2.00 x 10' 12 

0.1 

4.6 x lo 17 

800 

3 x lo' 5 

1.47 x 10' 13 

.01 

1.1 x 10 19 

l6o 

1 x 10 -5 

1.37 x 10 x 


TABLE II 

CALCULATED PERFORMANCE OF ADVANCED SILICON SOLAR CELLS 



Present 

Cells 

ADVANCED CELLS 

Bulk Resistivity, p, fl-cm 

10 

10 

l 

0.1 

.01 

2 

Short-Circuit Current, ^ sc > A/cm. 

.037 

.0442 

.0441 

.0433 

.0399 

Open Circuit Voltage, V QC > v 

.545 

.550 

.616 

.683 

.742 

Fill Factor, F.F. 

.72 

.792 

.806 

.817 

OO 

ro 

-3 

Efficiency, io 

10.7 

l4.i 

16.1 

17-8 

18.0 
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